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We developed the quasi-particle random phase approximation (QRPA) for the 
neutrino scattering off even-even nuclei via neutral current (NC) and charged cur- 
rent (CC). The QRPA has been successfully applied for the /3 and (3(3 decay of 
relevant nuclei. To describe neutrino scattering, general multipole transitions by 
weak interactions with a finite momentum transfer are calculated for NC and CC 
reaction with detailed formalism. Since we consider neutron-proton (np) pairing as 
well as neutron-neutron (nn) and proton-proton (pp) pairing correlations, the nn -|- 
pp QRPA and np QRPA are combined in a framework, which enables to describe 
both NC and CC reactions in a consistent way. Numerical results for u—^'^C, — ^^Fe 
and — ^^Ni reactions are shown to comply with other theoretical calculations and 
reproduce well available experimental data. 

PACS numbers: 21.60.Jz,23.40.Hc 

I. INTRODUCTION 

Neutrino (z/) (antineutrino (u)) scattering with a complex nucleusplays important roles 
on understanding the nuclear structure probed by weak interaction [l|-l3| as well as relevant 
u parameters in the u physics 4, sl, such as mass hierarchy and mixing angle 6'i3 through 
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detailed analysis of nuclear abundances in the core collapsing supernova (SN) explosion. 
Recently, lots of interests have been focused on the v process j4-7| in the nuclesynthesis 
because the emitted neutrino flux is expected to be large enough to excite the relevant 
nuclei in spite of the small cross section due to the weak interaction. Therefore cross 
sections for the neutrino(antineutrino)-nucleus [y{y^ — A) scattering with relevant nuclei 
are to be treated as important input data for a network calculation estimating the nuclear 



a 



abundances, in specific, for the weak rapid process 

Incident v{p^ energies exploited in the SN explosion [4, 5| are focused on the energy 
range from a few to tens of MeV because relevant z/(z/) energy spectra emitted from 
a proto-neutron star are presumed to be peaked mostly around the energy region by 
'oUowing a Fermi-Dirac distribution given by temperature (T) and chemical potential [a) 

1 El 

^^^"^ " F(a)T3 exv\{E,lT) - a] - 1 ' 
where Fia) is a nomalization factor. (T, a) can be chosen for neutrino types. Therefore, 
the v{y) — A reactions are usually sensitive on the collective motions of nucleons inside 
nuclei and proceed via two-step processes, i.e. the target nuclei are excited by the incident 



v{y) and decayed to the lower energy states with the emission of some particles {q]. 

The excitation occurs through various transitions, i.e. super allowed Fermi (J'^ = 0"*"), 
allowed Gamow Teller (GT) (J'^ = 1+), spin dipole (J'^ = 0^, 1^,2^), and other higher 
multipole transitions. Therefore dominant contributions of the two-step process stem 
from discrete and giant resonance (GR) states of the compound nucleus. Their typical 
excitation energy is below tens of MeV. It should be noticed that not only isobaric analogue 
state (IAS) and GT transitions but also contributions from higher order multipoles could 
play important roles in the v — A reaction p, 11Q|. 

In the experimental side, only a few data measured from accelerated-based v facilities, 
which are deduced mainly from ^^C, exist exclusively for the first stage of the two step 



processes, that is, the formation of compound nuc 



averaged total cross sections in the last decade 
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ei. They have been reported as v flux 



Hl5| 



JdE,a,{E,)f{E,) 
JdEJiE.) ' 

where a^{E^) and f{E^) are v — A cross section and corresponding neutrino flux. Unfor- 
tunately, there are no data including the decay processes, that is, the second stage in the 
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two-step process. Neutrino facilities usually used neutrino source from the pion decay at 
rest (DAR) and in flight (DIF). For the DAR neutrino, one obtains continuous spectrum 
for i^e; which is given as the Michel spectrum 

f{E^) = ^{m,-2E,) , (3) 



and different continuous spectrum for z/^ ll|, and mono-energetic u^^ fixed as 29.8 MeV. 



Continuous u^^ spectrum is obtained from the DIF neutrino. Recently, another new 
method for the production of intense and pure u beams has been proposed, called the 



beta-beam 
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17j ]. which is originally planned to study the CP violation in the lepton 
sector. The concept of the beta-beam is designed by using boosted radioactive ions which 
can be decayed with emitting neutrinos. For the study of fundamental interactions and 
nuclear structure, this new future facility may produce monochromatic u beams through 
the electron capture of radioactive ion beam. Of course, these beams depend on the exis- 
tence of ions decaying through the electron capture. It was also suggested how to obtain 



a mono-energetic u beam by using an electron capture process on a nucleus 18|, Il9 |. 

Most data are focused on ^^C target because the ^^C nucleus itself serves as a z/(z/) 
detector. Data for the inclusive reaction, ^^C(z/e, e~)^^N*, show about 3.3 ~ 6.7 in the 
10~^^ cm^ unit including indicated experimental error bars, while data for the exclusive 
reaction like ^^C(z/e, e")^2Ng.,. and ^'^C{u,u'Y^C (15.11 MeV) for the DAR neutrino are 
restricted to 6.07 ~ 10.4 and 8.5 ~ 12.3 in the same unit, respectively fll-[l5|. 



The u scattering data based on accelerator are mainly restricted to the formation of 
compound nuclei. But the excited states in compound nuclei subsequently decays into 
other nuclei with emitting particles such as proton, neutron, alpha, 7, and so on 
To describe the decay processes at the second stage, one needs additional calculations 
or branching ratios into the decay processes like Hauser-Feshbach (HF) statistical model 
3], [2o|. Any final state interactions (FSI) between outgoing particles and residual nuclei 
should be also taken into account at this stage. Of course, for the nuclear abundance, 
both two stages are to be successively considered. 

Before presenting our results, we briefly summarize recent theoretical status about 
u — A scattering with detailed references. Since the pioneering work by J. S. Cornnell et 
a/.Q, which exactly predicted the z/— ^^C reaction via char ged, current 20 years ago before 
the advent of the data, many theoretical calculations! 
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25| have been done 



and compared to the experimental data, mainly to the v —^"^ C reaction. Conventional 



approaches for the u — A scattering in the low energy region are shell modelfSM) jsf, 
random phase approximation (RPA) 2l|, and Quasi-particle RPA (QRPA) lOl. |21|. 

21| converge more or less on the experi- 



Alt hough results of most SM calculations 
mental data, they inevitably depend on the partic 
the so called qa quenching, and so on. The RPA 



e model space, the give n Hamiltonian, 



2l| and the QRPA (lol l2l| calculation 



usually overestimated the data 

[ 

culations are presented at Ref. [9 



jy a factor of 4 ~ 5. More detailed summaries of the cal- 



23| . But recent calculation by proton- neutron QRPA 



10| shows more improved results by simultaneously considering other relevant processes. 

In medium heavy or heavy nuclei, neutron-proton (np) pairing as well as proton-proton 
(pp) and neutron-neutron (nn) pairing could be important because of the small energy 
gap due to the neutron excess, although its empirical magnitude is small compared to 
those of light nuclei 2J]. Moreover the np pairing enables us to consistently describe the 



CC and NC reactions in a framework, as shown later on. 
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25 



includes continuum 



On the other hand, the Continuum RPA (CRPA) |22|, 
excitation spectrum in the quasi-elastic (QE) region in addition to the discrete and GR 
states, while most of SM, RPA and QRPA do not take explicitly the contribution from 
the QE region into account. The CRPA shows almost identical results compared to other 
approaches for the v — A scattering by the DAR neutrino, so that the QE contribution 
seems to be small enough to be neglected on a few tens of MeV region. 

In this work, we present our QRPA formalism for the v—A reaction via CC and NC type 
by including nn, pp and np pairing and compare them with the available experimental 
data and other calculations. Our QRPA has successfully described /3, 2z//3/3 and Oz/2/3 
decays 2J]. Fairly reasonable results with the experimental data are obtained. In Sec. II, 
we presented our whole formulas ioi v — A scattering, which is based on the Bardeen 
Cooper Schriffer (BCS) theory for the ground state description and the QRPA for the 
transition to the excited state. Relevant weak transitions and their operators are also 
summarized in Sec. II with the cross section formula. Numerical results and discussions 
are done in Sec. III. Summary and conclusions are given at Sec. IV. 
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II. FORMALISM 



A. Description of ground state with np pairing 



We start from the following hamiltonian 



H = Ho + H, 



int 



Ho 



S (^aa^'aa^o-ct i Hint S ^aabficydS^aa^blS^d.^^^'y i (4) 



abed ,a/5j5 

where the interaction matrix V is the antisymmetrized interaction with the Baranger 
hamiltonian 2^] in which two ^ factors, from J and T coupling, are included, so that 
the Hint in Eq.(4) is equivalent to the usual Hint = \T.a'b' c' d' ^ap-yS^aabpc-idsctaCtpCdsCc^- 
Roman letters indicate the quantum numbers of the nucleon states, i.e. a = {na,la,ja)- 
But Roman letters with prime include also the magnetic quantum number i.e. a = 
{ua, la,ja, '^a)- The isospiu of real particles is denoted as a Greek letter while the isospin 
of quasi-particles is expressed as a Greek letter with prime. The operator c'^^, Caa stand 
for the usual creation and destruction operators of the nucleons in the state of a with an 
isospin a. cj§ = Cn^./dj^.-m^ is the time reversed operator of Cds- 

We transform this hamiltonian by the general HFB transformation to quasi-particles. 



Cfl 



V(^c'^,'d'545 + ^c'^,'d'5Cd5) , a-c^,' = S «-'^' j'^Q-^ + t;*,^,^-,^c+ ) 



(5) 



d'S ~ ' ~ ' ~" d'5 

Here time reversal and spherical symmetries for the nuclei are presumed, so that we do not 
mix different single particle angular momentum states in Eq.(5). Then the Hamiltonian 
can be represented in terms of the quasi-particles as follows 



H =H,+ j:E^^>a+^,a^^> + H,^,,nt • (6) 

a a 

Using the transformation of Eq. (5) we obtain the following HFB equation 



6p \p 



A 



pp 
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(7) 





-^n ^np 

^pp ^pn f p ~l~ ^p 

\ Anp Ann (in ~\~ / 

, where E^^i is the energy of a quasi particle with the isospin quantum number a in the 
state c. If we neglect A„p, this equation reduces to the standard BCS equation. The 
pairing potentials Ap, A„ and Ap„ in Eq.(7) are detailed at Ref 24 1. 



B. QRPA coupled by neutron-proton pairing 



Excited states, |m; J'^M), in a compound nucleus, are generated from the ground state 
of the target nucleus, which is an even-even nucleus and acts as the BCS vacua for the 
quasi-particle, by operating the following one phonon operator to the initial nucleus 



Q 



JM ^ ,l^-(kf,'lu'j) 
klfj, u 



where 

C{kfi'lu'jM) 



S C Z^,,m,a,,'ak,' , C{kf, lu JM) = {-y-'''C{kf, lu J ~ M) (9) 



with a quasi-particle annihilation operator a^^i and Clebsh-Gordan coefficient Cj^^^^^^. 
If neutron-proton pairing is neglected, the phonon operator decouples to two phonon op- 
erators. One is for the charge changing reaction such as beta decay and CC neutrino 
reaction. The second is for the charge conserving such as electro-magnetic and NC neu- 
trino reactions. The amplitudes Xaa,bii and Yaa,biii which stand for forward and backward 
going amplitudes from state aa to 6/3, are obtained from the following QRPA equation, 
whose detailed derivation was also shown at Ref. 2^ 
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The A and B matrices are given by 



' ^{albl)J * 
vm 

^{a2b2)J 
vm 

^{alb2)J 



ym 

^ {albl)J 



{a2b2)J 



\ ^{alb2)J / 



' ^{klll)J * 
ym 

^(k2l2)J 
ym 

^(kll2)J 



•ym 

^ {klll)J 
ym 

^ {k2l2)J 

\ ym I 

\ -^(fcl«2)J / 

(10) 
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ap-yS 

where u and v coefficients related to the occupation probabiUties of a given state are 
determined from the HFB calculation in Eq.(7) with the pairing stren gth gnn, Qpp and gnp 
adjusted to the empirical pairing gaps A„„, App and 5„p, respectively [24] • E^^^i indicates 



the quasi-particle energy of the state k with the quasi-particle isospin 7', and (J^^'i^i = 1 
if A; = / and 7' = 5', otherwise cr^^';^' = a/2. The G{F) matrices are two body particle - 
particle (hole) matrix elements obtained as solutions of the following Bethe - Goldstone 
equation, called as Q matrix, 

y{W)ab,cd = Vab,cd + Kb,cd 7}-^/ (wjab,cd , [i-^) 

w - Ho 

where a, b, c, d indicate the single nucleon basis states characterized by oscillator type 
wave functions with single particle energies from the Woods-Saxon potential. Hq is the 
harmonic oscillator hamiltonian and Qp is the Pauli operator. V^^f^ is the phenomeno- 
logical nucleon-nucleon potential. We used the one boson exchange potential of the Bonn 
] . Of course, one can use effective interactions for the 2-body interactions, such as 



group 



the effective Skyrme force and the surface delta interactions [10| instead of the ^-matrix 



C. Description of CC and NC reaction 

Under the second quantization, matrix elements of any transition operator O between 
a ground state and an excited state \u; JM > can be factored as follows 

<QRPA\\dx\\uj;JM >=[Xy^j:<a\\dx\\b><QRPA\\[ctcb]^\\io;JM > . (14) 

ab 

Here the ffist factor < a| |& > can be calculated model independently for a given single 
particle basis 29|- Ground and excited states developed in the previous subsection are 
exploited for the second factor with the Quasi Boson Approximation (QBA). By using 
the phonon operator Q^'l^ in Eq.(8), we obtain the following expression for NC and CC 
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neutrino reactions. For NC reaction 



<QRPA\\dx\\uj;JM > (15) 
S ^ [Af^^' i,f^' < aa'\\dx\\bp' > [n^^^' v^^^' X^^> j,^' + v^^^' Up^^' Y^^> ] 

aa 6/3 

where the nomahzation factor is given as -^aabii'^J) — 

1 — ^ab^Q,'/?' (— + ^afe^Q,'/?')- Without the np pairing correlation, this ex- 
pression can be reduced to the following simple form 

<QRPA\\dx\\uj-JM > (16) 

= S [hfapbp < ap\\dx\\hp> [UpaVpbXapbp + VpaUpbYapbp] 
ah 

_(^_ya+jb+Jj^^^^ <:bp\\0^\\ap> [upbVpaXapbp + VpbUpaYapbp]] + {p^n) , 

where summations for a and are performed in Eq.(15) since quasi particle a and 
(= 1,2) include the np pairing. On the other hand, for CC reaction, they are given as 

<QRPA\\dx\\uj;JM> (17) 
S ^ [A/;^' b^' <aa\\dx\\b/3' > [u^^^, v^f^^^' X^^, + v^^^ u^^^^ Y^^ ] . 

aoL 6/3 

This form is also easily reduced to the results by pnQRPA without pn pairing 

< QRPA\\6x\ \ iO]JM >= S [Mapbn < ap||(9A||6n > [UpaVnbXapbn + VpaUnhyaphn] ■ (18) 

aphn 

Aa a simple application of the above form, we can calculate decay, whose operator is 
defined as 

Pi,, = < pmn > [c^Cr^Uf , /5|m = WlAif = H^'P^-M , (19) 
as follows in our QRPA formalism 2J], 



(l;t.l|rV||0+) = }X-^'b,'n^a.'pVb,'n + yPb,')l^a.'p^b,'n]C^^^ (20) 

aba 13 

(0+||r+a||l+ ) = S J^(T;i'W)i^a«'pVn + ^(r^''a.')i^a«'pVn]C'oe/(6,a) 

aba /3 

with 

Coef{a,b)= 2 {{L^)j,\\a\\ikl)jb) . (21) 



D. Neutral and Charged Current Operators 



For nuclear weak current operators, we start from a weak current on the nucleon level. 
The weak current operator takes a — A'^ current form by the standard electro- weak 



theory, which has isoscalar and isovector parts for the NC interaction 24 1 



= ^3^^ - - 2sm'ewJL - liVs' - A^) (22) 
= (1 - 2sm^ewM - ^3 - 2sm''ewV,^ - ^{Vj" - A^:) , 

with Weinberg angle 9w, where we used J^^ = + KT- ^3^(= ^1+12) "^si— ^1+12) 
are plus components of the isovector V^^ and by isospin rotation. Strangeness contri- 
butions, which are isoscalar parts, could be considered in —^{Vf — A^). 

For the CC interaction, only V^—A'^ term is involved, while = V3'^+V^'^ is concerned 
with meson electro-production. Therefore the CC reaction of the i^(z^) scattering is nearly 
independent of the strangeness contents. For the elastic scattering of polarized electron 
on the nucleon, J'' = — 2sin^6'vi/<^i'm — is exploited. 

For a free nucleon, the current operator comprises the vector, the axial vector and the 
pseudo scalar form factor, F^i^Q"^), Fa^Q"^) and Fp{Q'^) 

= F^iQ')r + ^2iQ')^^"'Q. + FAiQ'h'-l' + ^^^^"7' ' (23) 

where we take the scalar and the tensor form factor to be zero because of current con- 
servation and no existence of second class current, respectively. By the conservation of 
the vector current (CVC) hypothesis with the inclusion of the isoscalar strange quark 
contributions Ff(Q^), and the vector form factors for protons and neutrons F^' ^^^\Q'^) 



are expressed as 



32| 



Fr^^^-\Q') = (i -2sm'e^)F[^-\Q') - If-^'\Q') - ^i^^(g^) for NC (24) 
= (Ff (Q2) - F^Q')) for CC . 



The axial form factor is usually given by 33 1 



F^'^iQ') = l{T9A + 9%)/{l + QVMlf for NC (25) 



Fk'^iQn = -9A/il + Q'/Mir for CC , 

where qa and Ma are the axial coupling constant and the axial cut off mass, respectively. 
— (+) coming from the isospin dependence denotes the knocked-out proton (neutron). 



respectively. Since the energy region considered here is below the quasi-elastic region, the 
strangeness contributions are not taken into account in this report. Since we take + sign 
for Fa{Q'^) in Eq.(23), the axial form factor in Eq.(25) is just negative to the form factor 



elsewhere, for example, in Ref. 



32|. 



E. Relevant Operators for weak interactions 

In order to calculate weak interactions with nuclei, one resorts to the Hamiltonian, 
'H(x) = J dx.WJ^'^"''^{x.)Af^{qx.), contracted with weak hadron (nuclear) current oper- 
ator, VF,y"^"*^(x) and external weak field, A^lgx) = li,exp(—iq ■ x) with q = |q| and lepton 



current [2|. Here we follow the kinematics at Ref. 30|. Four momentum transfer is de- 



fined as g^(= qf^ = {k^ — u^)^) = q^ — c^ = —Q^ < 0, where and z/^ are final and initial 
lepton momenta. The excitation energy of a compound nucleus is given as w = —go > 0. 

The weak field can be expanded in terms of multipole operators by using two basic 
operators 

^(gx) = jj{q^)Yj'-^m , Mf/(gx) = j,(?x)Yf/,((].) , (26) 

where vectorial spherical harmonic Y^^l^iVLx) is expressed in term of spherical harmonic 
y£"((]^). I.e. yJ2i(^x-) = SmA < Lml\\{Ll)JMj > Y['{n^)ex with the spherical unit 
tensor e^. If we make non-relativistic reduction of the one body nucleon operator using 
Dirac wave function, we reexpress the one-body weak transition current in terms of 4 
different transition operators (Coulomb, longitudinal, electric and magnetic) as follows 

(T) 

Mjm;TmAi^) = [Ff )M^(gx) - ^±[F^nf'iq^) + ' i:" f M]]I^%?7) 

Cjm;TmA<1^) = [^Ff)M^(gx)+.(Ff)--|^Ff))S"f^(gx)]/--, 

q ZIVl]y 

ffM.,TMM^) = [-|[Ff)A'f^(gx) + l/.(-)S-^(gx)]+.Frs"/'(gx)]/-- , 
TZ^mM^) = -^[[^rAf^(gx)-i;.(-)E'f-'(gx)] + Ff)s-^(gx)]/--, 

where the superscript T(= 0, 1) means isoscalar and isovector. The 8 relevant single 
particle operators {Mf^^lf.^i'l^T.f.AfXf'Xf.^f') are given as follows. 
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whose detailed derivation was done at Ref. 



M 
J ' 



(2^ 



,M, 



1. 



MS7(gx).-V, 

[J]-^[-J^M*5V,(gx) + (J+i; 



jj-i(gx)] ■ -V 



spMj 



A/,, 



[J]-^[- J^My/_,,(gx) + (J + l)^/^Mji'_,(gx)] ■ a 

nf^(gx) = Mf^(gx)a.iv , fi'/^(gx) = ^^^(gx) + l/2S"/^(gx) . 

In actual calculation, the last operator is not used because it can be combined by other 
operators 10]. Single nucleon form factors in Eq.(27) are denoted as F^\q'^) with 
T = 0, 1 and X = 1,2,A,P standing for Dirac (X=l), Pauli (X=2), axial, and pseudo- 
scalar form factors, respectively. Detailed form factors are referred from Ref. {2!, 29| • The 
induced pseudoscalar form factor is usually parameterized by the Goldberger-Treimann 
relation 



2M, 



N 



(29) 



g2 + m2- 

where is the pion mass. The contribution of the pseudoscalar form factor vanishes for 
the NC reaction because of the negligible final lepton mass participating in this reaction. 
But it can contribute to the CC reaction, in particular, for z/^j — A reaction. 



F. Cross Sections 



Based on the initial and final nuclear states, cross section for ^{v) — A scattering 
through the relevant transition operator in Eq.(27) is given as 30 1 



,dai, 



Glek 



TT {2Ji + 1) L J=0 



[ Eil + u-P)\<Jf\\Mj\\J,> 



(30) 



+ (l-z7./3 + 2(z>-g)(g-/3))| < Jf\\Cj\\Ji > \ - q ■ {u + f3)2Re < Jf\\Cj\\Ji > < Jf\\Mj\\Ji > 

+ E (1 - (i> ■ m ■ m < jfWffu > r + 1 < jfiifT^'u > T) 



j=i 



± S g ■ (i> - /3)2i?e[< Jf\\7T'''\\J^ > < Jf\\rf\\Ji > 



j=i 



where (±) means the case of i^(z^), respectively, u and k are incident and final lepton 
3-momenta, and q = k — u, P = k/e with the final lepton's energy e. Of course, the 
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extremely relativistic limit (ERL) may yield more simple formula, but we use the general 
expression for the — A reaction. For CC reaction we multiplied Cabbibo angle cos^Qc 
and include the Coulomb distortion of the outgoing leptons due to the residual nucleus 

Since the Fermi function, FiZ ± l,e/), is deduced for the outgoing electron in the /3 
decay [s^, which usually assumed s-wave electron, one needs more deliberate approach 
for V — A reactions, in which outgoing lepton energy is higher rather than that of the /3 
decay. To exactly describe the Coulomb distortion, one needs to solve the Dirac equation 
of the outgoing lepton under the Coulomb potential due to the residual interactions with 
daughter nucleus (35|, l36|. But, since the exact solution needs a time consuming compu- 



tation, one usually exploits the averaged momentum of the outgoing particles. Typical 
method is the effective momentum approach (EMA), which is able to reproduce exact 



calculations more or less 



10 



37 



]. Here we use both approaches, the Fermi function 
and the EMA approach. By following the prescriptions on Ref . 10|, l21[ , below the energy 
on which both approaches predict same values, the Fermi function is used, while the EMA 



is adopted above the energy region. 



III. RESULTS 

A. v{p) - "^^C reaction 

In Fig.l, we show the cross sections of ^^C(t'e; e~)^^N* reaction for J'^ = 0^ ~ 4^ 
states with dominant transitions (0"'", 1^ and 2~ states' contributions) and total sum of 
all states. Other states' contributions not indicated here are within only a few percentage. 
Experimental LSND data for the exclusive ^^C(z/e, e~)^^Ng <..(i+) reaction are taken from 



Ref. HI]. Our results for the Gamow Teller 1"^ state clearly reproduce the data. Other J'^ 



states' contributions including spin dipole resonances turned out to become also important 



10|, in specific, beyond 45 MeV energy 



for understanding the v—^'^Q reaction via CC 
region. 

Results for the NC reaction, ^'^Q{y,v'f'^Q\ are shown in Fig. 2. Remarkable point is 
that the GT transition dominates overwhelmingly the cross section by NC. Below the 55 
MeV region, other transitions contribute within only a few %. This is contrast to those for 
CC reactions. Energy folded (or flux averaged) cross sections by Eq.(2) for the exclusive 
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and the inclusive reaction via CC and NC reactions are tabulated in table 1 with other 
theoretical results and experimental data. 

In order to pin down the ambiguities from the nuclear structure, results for /3 decay 
are also shown in table I. Here we used gA = 1-23 from Ref. which reproduces half 
lives of P~ and P~^{EC) decays, as shown in table 1. For pairing interactions, Qnn = 
1.6149, gpp = 1.4988 and gnp = 2.0698 are adjusted to fit the empirical pairing gaps 



Mi- 



Ann = 4.548, App = 4.430 and 6np = 2.489 MeV, respectively 

Our QRPA results reproduce most of the experimental data and show reasonable con- 
sistency with other theoretical model calculations. But the < cr > for the exclusive 
reaction, ^^C(z/e, e~)^^Ng,s., in table 1 overestimates the data about 10 % maximally. It 
might be understood from the following discussion about the Coulomb distortion. 

Main contributions to the folded cross section of the exclusive reaction come from the 
energy region E^^ < 52 MeV because of the Michel spectrum, Eq.(3). We divide the energy 
region into two regions for the Coulomb correction. In the energy region below 45 MeV, 
we take the Fermi function. But the Gamow Teller 1"*^ state transition, ^^C(i^e, e~)^^Ng.s,, 
overestimates the experimental data about 7 ~ 8 % as shown in Fig.l. This is the reason 
for the overestimation of the folded cross section. 

Here we discuss the Coulomb distortion effect. In Fig.3, we show three different results 
; Fermi function correction, effective momentum approach (EMA) [lO| and no Coulomb 



MA correction 
But, in the 



correction. The Fermi function correction turns out to be larger than the 
in the energy region above 60 Mev, which confirmed discussions in Refs. jol. 
energy region 4:5 < < 60 MeV, results by both approaches are nearly indiscernible and 
well reproduce the experimental data. However, below 45 MeV, the Fermi function and 
the EMA overestimate the cross sections about 7 ~ 8 % and 10 ~ 12 %, respectively, while 
the case of no Coulomb correction seems to match with the data. Therefore, for more exact 
result within 1 or 2 % for the flux averaged cross section for the ^^C(z/e, e~)^^Ng niore 
careful treatment of Coulomb corrections in the low energy region would be necessary. 

Our QRPA includes not only proton-proton and neutron-neutron pairing but also 
neutron-proton (np) pairing, but the contribution by the np pairing turns out to be 
only within 1 ~ 2 % for the relevant weak interaction in ^^C, such as (3^ decay and the 
u—^'^C reaction. Since the energy gap between neutron and proton space is relatively large 
in light nuclei, the neutron-proton pairing is expected to be small. But in the medium- 
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heavy or heavy nuclei, the effect could be larger as shown in the results for z/— ^^Ni(^^Fe) 
reactions. 



B. - ^^Ni and - ^^Fe reaction 

Here we calculate reactions of ^^Ni(z/e, z/g)^^Ni*, ^^Ni(i/e, i^l)^^Ni* via NC, and 
^^Fe(z/e, e~)^^Co* via CC. Results for z/(z/e)— ^^Ni by NC are shown in Fig.4 and 5. Re- 
markable difference between z/g— and Ue —^^ Ni is the magnitude of cross sections, in 
which cross sections by z/g reaction are about 40 % larger than those by reaction. But 
the cross section is dominated by the GT transition. Contributions by higher multipole 
transitions are shown to be within only a few %, similarly to the case of ^'^C{ue, u'^y^C*. 
Higher multipole contributions to the NC reaction turns out to be much smaller than 
those to the CC reaction. It is also confirmed in ^^Fe(z/e, e~)^^Co* reaction as shown in 
Fig.6. 

Fermi function for the Coulomb corrections is used on the energy region below 
45 MeV and the EMA is taken beyond 45 MeV. For pairing interactions, Qnn = 
1.3134(1.0244), ^pp = 1.1558(1.0340) and Qnp = 1.2785(1.3057) are adjusted to the em- 



pirical pairing gaps A„„ = 2.148(1.425 



for ^^Ni (^^Fe) nucleus, respectively 24 1. 



App = 2.078(1.572) and 6np = 1.102(0.336) MeV 



Since we have an experimental data for ^^Fe(z/e, e~)^^Co* by the DAR neutrino given as 
energy weighted cross section, which is the only experimental u — A reaction data beyond 
^^C, we compared our results to the experimental data. If the np pairing is switched off, we 
obtain < a > = 141.8 xlO~^^cm^. But the np pairing enhanced it to 173.5 xlO' 



■42 



cm 



which is located within the experimental data 256±108±43 xlO ^^cm^. As discussed 



in Ref. [10|, previous QRPA calculations without np pairing could not reproduce the 
experimental data. 

In Fig. 7, we show GT strength feasible from the (p,n) or (n,p) reaction on ^^Fe or ^^Ni 
target. They are calculated as 

BiGT^) = < /IISa^r.^lK > |' ■ (31) 

Results in Fig. 7 did not use the quenching factor. Our results for total GT strength with 
the quenching factor fq = 0.74 are 11.38 and 4.41 for B{GT\). They are consistent with 
those of experimental data, 9.9 ± 2.4 39| and 2.8 ± 0.3 j40|, respectively. As shown in 
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Fig. 7, contributions from the higher excited states lead to a bit larger value for B{GT+) 
than the data. Similarly to the results by the shell model jo], which is a hybrid model i.e. 
shell model for GT and IAS transition and QRPA model for high multipole transitions, 
our results reproduced the experimental data by using the same quenching factor as the 
shell model. We expect that more fruitful data for the relevant weak transitions, which 
enables to constraint the ambiguities from the nuclear structure. 

For higher multipoles we did not use the quenching factor by following the discussions 
at Ref. jo]. But larger contributions compared to those by Ref. {g] are obtained for the 
higher multiople transitions. In Ref. Ig], 80 % of the data is explained by the GT and IAS 
transition. Therefore our results explain systematically the u—^^Fe and u—^^Ni reaction 
as well as the GT strength in the single /3 decay in a systematically organized QRPA 
model. For the supernovae application, we show the cross sections averaged by the SN 
neutrino spectrum in Eq.(l) in Fig. 8. Cross section via CC are much larger, about 3 
times, than those by NC, irrespective of nucleus species. Application to the z/ process for 
other relevant nuclei are in progress. 



IV. SUMMARY AND CONCLUSION 



We applied the QRPA to the u — A reaction by including multipole transitions up to 
J'^ = 4^ with explicit momentum dependence and evaluated the //-reaction to ^^C, ^^Ni 
and ^^Fe target. Our results for the z/(z/) — A reaction show quite consistent results with 
available experimental data and other theoretical calculations. 

Remarkable points from our results are summarized as follows. Firstly, not only GT, 
IAS and spin dipole transitions but also other multipole transitions could play roles for 
V — A reaction via CC below the quasi elastic peak. But NC reaction is dominated by 
the GT transition. 

Secondly, np pairing as well as nn and pp pairing correlations act as important ingre- 
dients for describing the BCS ground state, and contribute to some extent to the cross 
section of the u — A reaction. In particular, it affects the results for medium heavy nuclei 
such as ^^Ni and ^^Fe, which was already expected from our previous results for 2z/2/3 and 
Oi/2/3 decays. 

Finally, for more exact evaluation within a few % for the flux averaged cross section 
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of the future would-be experimental data, more deliberately chosen Coulomb correction 
would be necessary. In particular, experimental data via NC would be desirable to pin 
down the ambiguities on the Coulomb corrections. 

The QRPA is a very efficient method to consider multi-particle and multi-hole in- 
teractions and their configuration mixing, which successfully described nuclear reactions 
sensitive on the nuclear structure, such as 2i/2/3 and Oi/2/3 decays. Therefore the am- 
biguities feasible due to the nuclear structure can be pinned down by reproducing the 
data related to the /3 and P(3 decay. The extension of our QRPA calculation for u — A 



reaction to explicitly include the deformation |4l| , which turns out to be so important for 
the exotic nuclei frequently appeared, is under progress. It enables us to perform various 
nuclear weak reactions for stable and unstable nuclei. 
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Table. I: Comparison of calculated and measured flux averaged cross sections for the z/— ^^C 
reaction in units of 10~^^ cm^, and half life time of neighboring nuclei. The cross sections are 
folded by the corresponding DAR neutrino spectra, where the Michel spectrum is used for V(. 
and energy is fixed at 29.8 MeV. "K" and "L" mean Karmen and LSND groups r_esults, 



respectively. Shell Model(SM) and Continuum RPA(CRPA) results are cited from Ref.|3|l and 
Ref.^], respectively. (9.834*) is a result with no Coulomb correction. 





i2C(i.e,e-)i2N,.,. 


i2C(i/e,e-)^2N* 


/?(-) : (12B(1+) ^i2c) 


Exp. 


8.9 ±0.3 ±0.9 [12J"L" 


4.3 ±0.4 ±0.6 [12l"L" 


23.6 ms 




9.1 ±0.5 ±0.8 [13J"K" 


5.1 ±0.6 ±0.5 []J"K" 




Ours 


11.53 (9.834*) 


6.1 


21.33 ms 


SM 


9.06 ~ 8.48 


5.22 ~ 4.87 




CRPA 


8.9 


5.4 








12C[(^„^;)]12C* 


P^+\EC) : (12n(1+) ^12(3) 


Exp. 


10.4 ± 1.0 ±0.9 [13]"K" 


3.2 ±0.5 ±0.4 [15]"K" 


11.0 ms 


Ours 


9.92 


3.60 


10.34 ms 


SM 


9.76 ~ 8.27 


2.68 ~ 2.26 




CRPA 


10.5 
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Fig. 1: (Color online) Cross sections of the ^'^C{i'(,,e )-'^^N* reaction for Jj^ = 0^ ~ 4^ states. 
Main multipole states (0^, 1^ and 2") contributions and total sum are presented. Cross section 



via J'^ = 1"*" state is compared to the LSND experimental data, C(fe,e ) N 



g.s.(l+) 



Hi. 



20 
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Fig. 2: (Color online) Cross sections of the "'^^C(z^e) z^e)^^C* reaction for = 0^ ~ 4^ states. 
Each multipole state contribution and total sum are presented. 
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Fig. 3: (Color online) Comparison of two different Coulomb corrections and no Coulomb cor- 
rection for the ^^C(z^e) e~)^^Ng J, (;^+-) reaction. Arbitrary normalized Michel spectrum for DAR 
neutrino is also presented to indicate the contribution by the incident neutrino energy bin to 
the flux averaged cross sections. Data point for the J'^ = 1"*" state is the same as Fig.l i.e. the 



LSND data for ^^C{ue,e~r^^g.s.{i+) 
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Fig. 4: (Color online) Cross sections of the ^^Ni(t'e, i/g)^^Ni* reaction for = 0^ ~ 4^ states. 
Each multipole state contribution and total sum are presented. 
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Fig. 5: (Color online) Cross sections of the ^^Ni(i/e, i^e)^^^i* reaction for J^- = 0^ ~ 4^ states. 
Each multipole state contribution and total sum are presented. 
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Fig. 6: (Color online) Cross sections of the ^^Fe(fe,e )^^Co* reaction for = 0^ ~ 4^ states. 
Each multipole state contribution and total sum are presented. 
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Fig. 7: (Color online) Gamow- Teller strength B{GT±) for ^"Fe 
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Fig. 8: (Color online) Temperature dependence of the energy weighted cross section, Eq. [2] for 
v—^'^C, -^^Ni and -^^Fe reactions, where neutrino spectrum for the SN, Eq.(l), is exploited. 
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